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Abstract 

Certain one-dimensional Fermi systems have an energy gap in the 
bulk spectrum while boundary states are described by one Majorana 
operator per boundary point. A finite system of length L possesses two 
ground states with an energy difference proportional to exp(— L/Iq) 
and different fermionic parities. Such systems can be used as qubits 
since they are intrinsically immune to decoherence. The property of a 
system to have boundary Majorana fermions is expressed as a condi- 
tion on the bulk electron spectrum. The condition is satisfied in the 
presence of an arbitrary small energy gap induced by proximity of a 
3-dimensional p-wave superconductor, provided that the normal spec- 
trum has an odd number of Fermi points in each half of the Brillouin 
zone (each spin component counts separately). 



Introduction 

Implementing a full scale quantum computer is a major challenge to mod- 
ern physics and engineering. Theoretically, this goal should be achievable 
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due to the possibility of fault-tolerant quantum computation |]T|. Unlim- 
ited quantum computation is possible if errors in the implementation of each 
gate are below certain threshold £|, |JJ. Unfortunately, for conventional 
fault-tolerance schemes the threshold appears to be about 10~ 4 , which is 
beyond the reach of current technologies. It has been also suggested that 
fault-tolerance can be achieved at the physical level (instead of using quan- 
tum error- correcting codes). The first proposal of these kind [fj was based 
on non-Abelian anyons in two-dimensional systems. The relation between 
quantum computation, anyons and topological quantum field theories was 
independently discussed in [0. A mathematical result about universal quan- 
tum computation with certain type of anyons has been obtained recently [§], 
but, generally, this approach is still undeveloped. In these paper we describe 
another (theoretically, much simpler) way to construct decoherence-protected 
degrees of freedom in one- dimensional systems ("quantum wires"). Although 
it does not automatically provide fault-tolerance for quantum gates, it should 
allow, when implemented, to build a reliable quantum memory. 

The reason why quantum states are so fragile is that they are sensitive 
to errors of two kinds. A classical error, represented by an operator crj, flips 
the j-th qubit changing |0) to |1) and vice versa. A phase error cr| changes 
the sign of all states with the j-th qubit equal to 1 (i. e. j-th spin down, if 
the qubits are spins) relative to the states with the j-th qubit equal to 0. 
It is generally easy to get rid of one type of errors, but not both. However, 
the following method of eliminating the classical errors is worth considering. 
Let each qubit be a site that can be either empty or occupied by an electron 
(with spin up, say, the other spin direction being forbidden). Let us denote 
the empty and the occupied states by |0) and |1), respectively. (Such sites are 
not exactly qubits because electrons are fermions, but they can be also used 
for quantum computation ||). Now single classical errors become impossible 
because the electric charge is conserved. Even in superconducting systems, 
the fermionic parity (i.e. the electric charge (mod 2)) is conserved. Two 
classical errors can still happen at two sites simultaneously, but this would 
require that an electron jumps from one site to the other. Such jumps can be 
avoided by placing the "fermionic sites" far apart from each other, provided 
the medium between them has an energy gap in the excitation spectrum. 

Obviously, this method does not protect from phase errors which are now 
described by the operators aja,-. To the contrary, different electron configu- 
rations will have different energies and thus will pick up different phases over 
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time. Even without actual inelastic processes, this will produce the same 
effect as decoherence. However, a simple mathematical observation suggests 
that the situation could be improved. Each fermionic site is described by a 
pair of annihilation and creation operators a,j, aj. One can formally define 
Majorana operators 

c 2j -i = a s + at, c 2j = aj Qj (j = 1, . . . , N) (1) 

which satisfy the relations 

c m = c m , cic m + c m ci = 25 lm (I, m = 1, . . . , 2N). (2) 

If the operators c 2 j-\ and c 2 j belonged to different sites then the phase error 
a ] a i = \Q + * c 2j-i c 2j) would be unlikely to occur. Indeed, it would re- 
quire interaction between the two "Majorana sites" which could be possibly 
avoided. Note that a single Majorana operator c 2 j_\ or Cj can not appear 
as a term in any reasonable Hamiltonian because it does not preserve the 
fermionic parity. Thus an isolated Majorana site (usually called a Majorana 
fermion) is immune to any kind of error! 

Unfortunately, Majorana fermions are not readily available in solid state 
systems. The goal of this paper is to construct Hamiltonians which would 
give rise to Majorana fermions as effective low-energy degrees of freedom. 
Surprisingly, this can be done even with non-interacting electrons. (Some 
interaction is actually needed to create superconductivity, but it can be ef- 
fectively described by terms like AajOk). The general idea is quite simple. 
An arbitrary quadratic Hamiltonian can be written in the form 

H = l - Ai m cic m {A* m = A lm = -A ml ). (3) 

l,m 

Its ground state can be described as "pairing" of Majorana operators: nor- 
mal mode creation and annihilation operators a^, a m , which are certain lin- 
ear combinations of q, come in pairs. (In this sense, an insulator and a 
superconductor represent different types of pairing). In some cases, most 
Majorana operators are paired up with an energy gap while few ones (lo- 
calized at the boundary or defects) remain "free". For example, unpaired 
Majorana fermions exist on vortices in chiral 2-dimensional p-wave super- 
conductors [ lOl, We will show that Majorana fermions can also occur at 



the ends of quantum wires. 
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1 A toy model and the qualitative picture 



We are going to describe a simple but rather unrealistic model which exhibits 
unpaired Majorana fermions. It attempts to catch two important properties 
which seem necessary for the phenomenon to occur. Firstly, the U(l) sym- 
metry dj i — ► e cij, corresponding to the electric charge conservation, must 
be broken down to a Z 2 symmetry, dj i— > — dj. Indeed, if a single Majorana 
operator can be localized, symmetry transformation should not mix it with 
other operators. So we should consider superconductive systems. The par- 
ticular mechanism of superconductivity is not important; we may just think 
that our quantum wire lies on the surface of 3-dimensional superconductor 
(see fig. [ip. The second property is less obvious and will be fully explained 
in Sec. ||. Roughly speaking, the electron spectrum must strongly depend on 
the spin. Here we will simply assume that only one spin component (say, f) 
is present. Q 




Figure 1: A piece of "quantum wire" on the surface of 3-dimensional super- 
conductor. 

Consider a chain consisting of L ^> 1 sites. Each site can be either empty 
or occupied by an electron (with a fixed spin direction). The Hamiltonian is 

Hi = y^f— w(a\aj + i + Oj +1 a,j) — /i(a]a.,- — |) + Aajdj+i + A*aj +1 ajY (4) 
j 

Here if is a hopping amplitude, fi a chemical potential, and A = | A|e i6> the 
induced superconducting gap. It is convenient to hide the dependence on the 
phase parameter 6 into the definition of Majorana operators: 

C2j-i = e % ?a,j + e" l 5aj, c 2j = -ie^aj + ie~^a] (j = 1, . . . , L). (5) 

1 It appears that only a triplet (p-wave) superconductivity in the 3-dimensional sub- 
strate can effectively induce the desired pairing between electrons with the same spin 
direction — at least, this is true in the absence of spin-orbit interaction. 
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In terms of this operators, the Hamiltonian becomes 

Hi = ^X](-Aic 2 j_ic 2 j + (w + \A\)c 2j c 2j+l + (-w + |A|)c2j-iC2j+ 2 )- (6) 

3 

Let us start with two special cases. 

a) The trivial case: \A\ = w = 0, fi < 0. Then H 1 = —fiJ2j{ a ] a j ~ \) — 
|(— fj) J2j C2j-iC2j- The Majorana operators c 2 j-i,c 2 j from the same 
site j are paired together to form a ground state with the occupation 
number 0. 

b) |A| = w > 0, jj, — 0. In this case 

#1 = iw^CjyC^+i. (7) 

j 

Now the Majorana operators 02^,02^+1 /rom different sites are paired 
together (see fig. 0). One can define new annihilation and creation 
operators dj = \{c 2 j +ic 2 j+i), a] = \[c 2 j —ic 2 j+i) which span the sites 
j and j + 1. The Hamiltonian becomes 2wY^j=i{0'}o'j — \)- Ground 
states satisfy the condition dj\ip) — for j — 1, . . . , L — 1. There are 
two orthogonal states |^o) and \tpi) with this property. Indeed, the 
Majorana operators V = C\ and b" = c 2 l remain unpaired (i. e. do not 
enter the Hamiltonian), so we can write 

- ib'b"\^) = |Vo>, -ib'b"\^) = (8) 



Cl C2 C3 C4 C2L-1C2L Cl C2 C3 C4 C2L-1 C2L 




a) b) 

Figure 2: Two types of pairing. 

Note that the state l^o) has an even fermionic parity (i. e. it is a superposition 
of states with even number of electrons) while has an odd parity. The 
parity is measured by the operator 

p=iK-ic 2j „ lC2j ). (9) 
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These two cases represent two phases, or universality classes which exist in 
the model. A subtle point is that both phases have the same bulk properties. 
In fact, one phase can be transformed to the other (and vice versa) by mere 
permutation of Majorana operators, 

Cm !-> C m+1 . (10) 

Such a local transformation (operator algebra automorphism) is usually con- 
sidered as "equivalence" in the study of lattice models. Yet the boundary 
properties of the two phases are clearly different: only the phase (b) has 
unpaired Majorana fermions at the ends of the chain. This is due to the fact 
that the operators Qy-i, C2j belong to one physical site while C2j, C2j+i do not. 
We may put it this way: one can not cut a physical site into two halves; if 
one could, both types of boundary states would be possible in both phases. 

Also note that the transformation (|10|) can not be performed in a con- 
tinuous fashion, starting from the identity transformation. From the math- 
ematical perspective, it means that one should have different definitions for 
"weak" and "strong" equivalence of lattice models . We will not touch such 
abstract matters here. 

Now we want to study the model at arbitrary values of w, \i and A. Let 
us begin with some generalities. Let N be the total number of fermionic 
sites in the system, for now N = L. The Hamiltonian has the general 
form (^). Hence it can be reduced to a canonical form 

i N N 

^canonical = ~ e mKnKi = 12 e ( 5 L 5 m ~ §) (e m > 0). (11) 

m=l m=l 

Here b' m ,b m are real linear combinations of c^j-x^c^j with the same com- 
mutation relations whereas a m = \{b' m + ib m ), a^ m = \ib' m — ib''^). More 

2 Nonlocal transformations can change the physical properties of the model even more 
dramatically. The Jordan- Wigner transformation C2j-i i— > crj 111= i u k-> c 2j t— > &j ITfc=i a k 
takes our model to a spin chain with xx and yy interactions and a z-directed magnetic 
field. Unlike ©, the Jordan- Wi gner transformation is well defined at the ends of the 
chain. However, this mathematical procedure falls apart in the physical context, as far as 
perturbations are involved. Indeed, the phase (b) has now an order parameter (a x ) =/= 0. 
External fields will interact with the order parameter breaking the phase coherence between 
\ip Q ) and |V>i). 
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where W is a 2iV x 2iV real orthogonal matrix (W 7 ^ = iy 5 
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(12) 
7) whose 
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rows are eigenvectors of A. The numbers e m > are one-particle excitation 
energies. However, it is more convenient to deal with a "double spectrum" 
{e m , —e m } since the matrix A has eigenvalues ±ie m . 

The bulk spectrum (energy vs. momentum) is given by 



e(q) = ±v/(2u> cosg + fi) 2 + 4| A| 2 sin q, 



— TT < q < TT. 



(13) 



We may conjecture that the phases (a) and (b) extend to connected domains 
in the parameter space where the spectrum has a gap. The signs of \x and w 
seem not to be important, so we actually expect that the phase (a) occurs at 
2\w\ < \fx\ while the phase (b) occupies the domain 2\w\ > |/i|, A^O. (The 
phase boundary is given by the equation 2\w\ = while A = 0, 2\w\ > \fi\ 
is a line of normal metal phase inside the domain (b)). 

To verify the conjecture, we need to find boundary modes. They corre- 
spond to eigenvectors of A localized near the ends of the chain. Due to the 
spectrum symmetry e 1— ► — e, zero eigenvalues can occur in a general position. 
If exist, such zero modes should have the form 



b' = Ej(a+a ; + + a'_xi,) c 2 j-i 



x± 



fi±^//j, 2 -Aw 2 + A\A\ 
2(w + \A\) 



(14) 

We will consider two cases corresponding to the expected existence domains 
of the two phases. 



a) If 2\w\ < \fj,\, we have 



> 1, \x-\ < 1 or 



[or a + , a 



< 1, > 1. 
" N can be non- 



Therefore, only one of the coefficients a' + , a 
zero, depending on whether the mode is to be localized at the left or at 
the right end of the chain. This makes it impossible to satisfy boundary 
conditions. So the supposed zero modes ( |I4T ) do not exist. 
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b) If 2w > \fi\, A ^ 0, we find that \x+\, < 1. Hence b' is localized 
near j = whereas b" is localized near j = L. There are also boundary 
conditions a' + + a'_ = 0, a'+x'+ L+1 ^ + a'!_xZ = 0, but they can 
be satisfied too. The zero modes b', b" are actually the same as the 
unpaired Majorana fermions discussed before. If — 2w > \fx\, A ^ 
then b' and b" change places. Thus the unpaired Majorana fermions 
exist in the whole expected domain of the phase (b). 

The above analysis is exact in the limit L — > oo. If the chain length L 
is finite, there is a weak interaction between b 1 and b" . (For definiteness, we 
will always assume that b' is at the left end of the chain whereas b" is at the 
right end). This interaction is described by an effective Hamiltonian 

H cS = l -tb'b", toce- L / l °, (15) 



where l 1 is the smallest of 



In \x A 



In \x. 



[note that both logarithms 



and 

have the same sign). Thus the energies of the ground states j^o) an d l^i) 
(see eq. (§)) differ by t. Note that it is not obvious any more which state 
of the two is even and which is odd. In the case — 2w > \fi\, the parity 
is proportional to (— 1) L . (This factor is the parity of the bulk part of the 
chain) . 

The effective Hamiltonian flT5|) still holds if we include small electron- 
electron interaction (a four-fermion term) into @). Indeed, the physical 
meaning of t is an amplitude for a fermionic quasiparticle to tunnel across the 
chain. In a long chain, this amplitude vanishes as e~ L ^ l ° if the bulk spectrum 
has a gap. 

Finally, we will discuss a role of the phase parameter 9 (A = e* e |A|). 
According to eq. (Bl), the Majorana operators c 2 j_i,c 2 j are multiplied by 
— 1 when 9 changes by 2-k. The physical parameter A is the same at 9 and 
9+2n, of course, but the ground states should undergo certain transformation 
as 9 changes to 9 + 2n adiabatically. Note that the transformation c m i— > 
— c m also occurs if one conjugates c m by the parity operator P. Within the 
effective Hamiltonian approach, P is the same as s(L)(—ib'b") (s(L) = ±1). 
Hence the adiabatic change of the superconducting phase by 2n results in 
the unitary transformation 

V = s(L)(-ib'b") : 7|Vo) = |V'o>, Wi> = -|V>i>- (16) 
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This is equivalent to transfer of an electron between the ends of the chain. 
Some physical consequences of this result will be mentioned in Sec. §. 



2 A general condition for Majorana fermions 

Let us consider a general translationally invariant one-dimensional Hamilto- 
nian with short-range interactions. It has been mentioned that the neces- 
sary conditions for unpaired Majorana fermions are superconductivity and a 
gap in the bulk excitation spectrum. The latter is equivalent to the quasi- 
particle tunneling amplitude vanishing as e~ L ^ l °. Besides that, it is clear 
that there should be some parity condition. Indeed, Majorana fermions at 
the ends of parallel weakly interacting chains may pair up and cancel each 
other (i. e the ground state will be non-degenerate). So, provided the en- 
ergy gap, each one-dimensional Hamiltonian H is characterized by a "Ma- 
jorana number" M. = Ai(H) = ±1: the existence of unpaired Majorana 
fermions is indicated as M. = — 1. The Majorana number should satisfy 
M(H' © H") = M(H')M(H"), where © means taking two non- interacting 
chains. 

Remarkably, the Majorana number reveals itself even if the chain is closed 
into a loop. This is handy as it eliminates the need to study boundary modes. 
Let H(L) be the Hamiltonian of a closed chain of length L ^> Iq. (H itself is 
a template which is used to generate H(L) for any L). We claim that 



P(H(L 1 + L 2 )) = M(H) P{H(Li)) P(H(L 2 )), 



(17) 



where P(X) denotes the ground state parity of a Hamiltonian X (assuming 
that the ground state is unique). 



Li 



L 2 

m>o-o-o-o j 



Li 

M>0-0-CK>0- 



h" 
"2 



L 2 

O-O-OQO 1 



b', 



b", 



a) b) 
Figure 3: Reconnecting closed chains. 



The following argument justifies eq. (|17D. An open chain of length L 
can be described by an effective Hamiltonian which only includes boundary 
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modes. If Ai(H) = — 1, there are Majorana operators b', b" associated with 
the ends of the chain. The parity operator P (see eq. ([|)) can be replaced 
by s(L)( K —ib'b"), where s(L) = ±1. Thus the fermionic parity of \ip a ) is 
s(L) (—1)°, a = 0, 1. If we close the chain, the effective Hamiltonian is 
H eS (L) = '^ub"b' . (We have chosen to write b"b' in this order because b" 
precedes b' in the left-to-right order on the loop, where they are next to each 
other). The parameter u represents direct interaction between the chain ends 
(unlike t from eq. (|15D), so u does not depend on L. The ground state of the 
closed chain is if u > and l^o) if n < 0. Hence 

P(H(L)) = -s(L) sgnu. 

Now let us take two chains, one of length L±, the other of length L 2 . 
There are two ways to close them up, see fig. [3|. Both cases can be described 
by effective Hamiltonians: 

Hrtil^QHrtiLz) = l -u(b'(b' 1+ b' 2 %), H^L.+L,) = % -u {b'[b' 2 +b'ib' x ) . 
It follows that 

P(H(Li)) P(H(L 2 )) = s(Li) s(L 2 ), + L 2 )) = - s(Lx) s(L 2 ). 

So the equation ( |I7D holds for M. = — 1. It also obviously holds for M. = 1 
because in this case there are no boundary modes to worry about. 

Computing the Majorana number in general (especially for strongly cor- 
related systems) may be a difficult task. However, the computation can be 
carried through for any system of non-interacting electrons. Consider a peri- 
odic chain of L unit cells with n fermionic sites (i. e. In Majorana operators) 
per cell, which totals to iV = nL fermionic sites. We will index the Majorana 
operators as q q , where I — 1, . . . , L, a — 1, ... , In. The Hamiltonian is 

H = jEE B ^( m ~ C *" C m/3 {BolpUT = B a p(j) = -B/3 a (-j)y 

l,m ot,f3 

(18) 

We assume that the chain forms a loop, so m — I should be taken (mod L). 

Eq. (0) is a special case of (^), so we will first find P(H) for the general 
quadratic Hamiltonian (H), assuming that the matrix A is not degenerate. 
The canonical form of this Hamiltonian ([11]) has an even ground state |0). 
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The transformation ( |P2"D can be represented as conjugation by the parity- 
preserving unitary operator U = exp(| AmQCm) if W has the form 



W 



exp(D) for some real skew-symmetric matrix D, i.e. if det W = 1. 



Otherwise, the transformation (12) changes the parity. Hence 

P(H) = sgndetW = sgnPf A. 



(19) 



We remind the reader that the Pfaffian Pf is a function of a skew- 
symmetric matrix such that (Pf A) 2 = det A. It is defined as follows 



Pf A 



2 N N\ 



sgn(r)A 



r(l),r(2) 



A 



t(2N-1),t{2N)- 



(20) 



t£S 2 n 



(Here S2N is the set of permutations on 2N elements). For example, 



Pf 



V 





-ai2 
-ai3 
-ai4 
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«13 Ol4 ^ 

«23 a 24 

a 34 

-a 34 J 



a 120'3A + 014^23 — Ol3<3.24- 



In eq. ( p!9|) we have used this property of the Pfaffian: 

~Pi{W AW T ) = Pi(A) det(W). 



(21) 



Now we are to compute the Pfaffian of the matrix B from eq. flTjjl) . First, 
we use the Fourier transform, 

B al3 (q) = J2e iq3 B a/3 (j), q = 2v3 (mod 2vr), fe = 0, . . . , N - 1. (22) 



The matrix 5(g) has these symmetries: 

B\q) = -B(q) = B T ( 



(23) 



The spectrum e(g) is a continuous real 2n-valued function on a circle (real 
numbers (mod 2tt)) given by the eigenvalues of iB(q). It has the symmetry 
e(— g) = — e(g). The energy gap assumption implies that e(g) never passes 
0. It follows that there are n positive and n negative eigenvalues for any q. 
Indeed, this is the case for q = due to the e h -e symmetry, hence it is 
true for any q by continuity. 
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It follows from eqs. (|22| ) and ( f2~T| ) that 

PfB = ( J] Pffl(g)) ( II detS(g) ] . (24) 

Remember that g is considered (mod 27r), so q = —q when g = or q = n. 
In the q ^ —q case, each {q, —q} pair is counted once. Note that det B(q) 
is a positive number since iB(q) has n positive and n negative eigenvalues. 
Hence 



sgnPfS = J] sgn(Pf5(g)) 



sgn(Pf B(0)) sgn(Pf B(n)) if L is even, 

sgn(Pfi?(0)) if Lis odd. 

(25) 

Finally, we get 



M(H) = sgn(Pf B(0)) sgn(Pf 5(tt)). (26) 

This very general equation can be simplified if superconductivity is a 
weak effect, i.e. |A| l e (0)|, |e(7r)|. Indeed, the right hand side of (f26|) 
makes perfect sense for a L r (l)-symmetric Hamiltonian 

Ho = \j2J:C a0 (m-l)ala mP (C af3 (j)* = C Pa (-j)), (27) 

l,m a, (3 

where a,/3 = l,...n refer to fermionic sites. The eigenvalues of C(q) (the 
Fourier transform of C) form a "single spectrum" €o(q)- The "double spec- 
trum" defined above is e(q) = ±e (g). It is easy to show that PfB(q) = 
det C(q) for q = 0, tt. Hence 

M(H ) = (-l)"W-"(o) j (28) 

where v(q) is the number of negative eigenvalues of C(q). Note that u(tt) — 
z/(0) equals (mod 2) the number of Fermi points on the interval [0,tt]. (A 
Fermi point is a point where eo(g) passes 0). In the most interesting case 
v(n) — i/(0) = 1 (mod 2), the Hamiltonian H has a gapless spectrum. So 
eq. (^) is only relevant in the presence of superconductivity, i. e. a small 
symmetry-breaking perturbation which opens an energy gap. 
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3 Speculations about physical realization 



Physical realization of an Ai = — 1 quantum wire is a difficult task because 
electron spectra are usually degenerate with respect to spin, so z/(0) and 
v(tt) are even. The degeneracy at q = and q = it can be lifted only 
if the time reversal symmetry is broken. Thus spin-orbit interaction does 
not help. External magnetic field could help, but the Zeeman energy gfis'H 
is usually small compared to other spectrum parameters, so z/(0) and //(tt) 
do not change. The situation may be different for charge and spin density 
waves which add fine features to the electron spectrum. Charge density waves 
(CDW) tend to occur at the wave vector = 2q F so that a gap opens at the 
Fermi level. In the presence of magnetic field, qp is slightly different for the 
t and I spin components, so it is possible that matches only one of them. 
The resulting spectrum is shown in fig. ^ in the g*/(27r) units. This scenario 
can be realized if |A| < Ecbw S gi^s'H- 

Another speculative possibility is to use midgap states at the edge of a 
two-dimensional p-wave superconductor [[12 . 



A e (q) 



7i q 



spin up 
spin down 



Figure 4: An electron spectrum in the presence of magnetic field and CDW. 

A quantum wire bridge between two superconducting leads (see fig. |a) 
could be used as an experimental test for Majorana fermions. When the phase 
parameter 82 in the right piece of superconductor changes by 2ir (relative to 
6*1), a fermionic quasiparticle is effectively transported to the junction region. 
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At the same time, the Majorana fermions at the ends of the wire switch from 
l^o) to j^i) or vice versa. If the quasiparticle stays localized, the junction 
parameters change. They change back when 9 2 changes by another 2tt. Thus 
the Josephson current is 47r-periodic as a function of 9 = 9 2 — 8\. In fact, 
it is more accurate to say that the Josephson energy Ej is 27r-periodic but 
2- valued, as shown in fig. ||b. The two levels may not quite cross at 6 = n due 
to a non- vanishing tunneling amplitude t oc e _L/// °, where L is the distance 
between the junction and the closest end of the wire. 



Ej(9) 




a) b) 
Figure 5: A Josephson junction made of quantum wire. 

Interesting phenomena can also take place in the simple layout shown 
in fig. [I]. Suppose that the superconducting island supporting the quantum 
wire is connected to a larger piece of superconductor through an ordinary 
Josephson junction. If the Coulomb energy is comparable to the Josephson 
energy, spontaneous phase slips can occur. Each 2n phase slip is accompanied 
by the operator V (see eq. [T^). The phase slips occur by tunneling, so the 
effective Hamiltonian is 

H eSA = -\V -\*V^ = -s(L)tb'b", t = 4ReA, (29) 

where A is the amplitude of the 9 i— ► 9 + 2ir process while A* corresponds 
to the reverse process. Similarly, if the superconducting island supports two 
quantum wires, the effective Hamiltonian becomes 

#eff. 2 = -AViV 2 - \*VM = ^(iiM^a) * W 2 &2- (30) 
Turning A on and off can be possibly used for quantum gates implementation. 
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